Can our Universe be inhomogeneous on large sub-horizon 

scales? 



P. P. Avelino 1 - 2 * J. P. M. de Carvalho 1 ^ C. J. A. P. Martins 4 ' 1 * and J. C. R. E. 

Oliveira 1 

1 Centro de Astrofisica, Universidade do Porto 
Rua das Estrelas s/n, 4150-762 Porto, Portugal 
2 Dep. de Fisica da Faculdade de Ciencias da Univ. do Porto 
Rua do Campo Alegre 687 - 4169-007 Porto - Portugal 
3 Dep. de Matemdtica Aplicada da Faculdade de Ciencias da Univ. do Porto 
Rua das Taipas, 135 - 4050 Porto - Portugal 
4 Department of Applied Mathematics and Theoretical Physics 
Centre for Mathematical Sciences, University of Cambridge 
Wilberforce Road, Cambridge CBS OWA, U.K. 
(February 5, 2008) 

We show that our Universe may be inhomogeneous on large sub- 
horizon scales without us being able to realise it. We assume that a 
network of domain walls permeates the universe dividing it in domains 
with slightly different vacuum energy densities. We require that the 
energy scale of the phase transition which produced the domain walls 
is sufficiently low so that the walls have a negligible effect on structure 
formation. Nevertheless, the different vacuum densities of different 
domains will lead to different values of the cosmological parameters 
QP A , and h, in each patch thus affecting the growth of cosmological 
perturbations at recent times. Hence, if our local patch of the universe 
(with uniform vacuum density) is big enough - which is likely to happen 
given that we should have on average about one domain per horizon 
volume - we might not notice these large-scale inhomogeneities. This 
happens because in order to see a patch with a different vacuum density 
one may have to look back at a time when the universe was still very 
homogeneous. 
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I. INTRODUCTION 



The last year or so has seen the first serious attempts to provide some direct 
connections between "fundamental" high-energy physics [Q and "low-energy" stan- 
dard cosmology Although this "top-down" approach is still at a very early stage, 
a number of crucial general trends already became apparent. For example, since 
high-energy theories are formulated in higher dimensions, any low energy limit will 
necessarily involve dimensional reduction, and possibly also compactification H. 
This turns out to be crucial because as a result of this process the low energy, 
four-dimensional coupling constants become functions of the radii of extra dimen- 
sions, which are often variables. One can therefore end up with low-energy effective 
models in which some of the "fundamental" constants of nature are time and/or 
space-varying quantities. There are a number of known examples of such models 
. On the other hand, there are recent tentative suggestions of a time- variation 
of the fine structure constant |J , but these require further confirmation. It is there- 
fore interesting to study the possible observational signatures of such variations, 
and in particular to find out how such observational signatures constrain the pos- 
sible models. It turns out to be easier to study this issue by constructing simple 
"toy models" . This provides a "bottom- up" approach, in which one gives up the 
possibility of testing particular assumptions from first principles, but instead has 
the possibility of exploring a larger patch of parameter space. This idea goes back 
at least to Dirac, and had its first detailed realization with the Brans-Dicke model 
JTo|] , which has a varying G. A number of toy models have recently been constructed 
to analyse possible variations of the fine structure constant [|ll-13|, the speed of 
light Jl|-H and electric charge [§§]. 

A somewhat related approach is that of "quintessence" (see for example p^-[27|]). 
These are essentially models with a time- varying cosmological constant. 

Here we consider the possible observational effects of having a universe made up 
of different domains, each with a different value of the cosmological constant. Such 
a structure would dramatically influence the future evolution of the universe J2H.29]. 
Recent observations of Type la Supernovae up to redshifts of about z ~ 1 |3(1-32], 
when combined with CMBR anisotropy data, seem to indicate that our patch of 
the universe is currently characterised by the parameters ~ 0.7 and £l® n ~ 0.3, 
implying that the cosmological constant become important only very recently. As 
has been pointed out before, the Supernovae measurements are local, and so they 
can not be extrapolated all the way to the horizon. For example, we could be living 
in a small, sub-critical bubble, and our local neighbourhood could have a value of 
~ 0.7 that is uncharacteristic. Here we discuss some basic consequences of such 
a scenario. We shall assume that different regions of space have different values of 
the vacuum energy density, separated by domain walls. This can be achieved if there 
is a scalar field, say (j>, which within each region sits in one of a number of possible 
minima of a time-independent potential. The above simplifying assumptions could 
be relaxed; for example one could instead consider quintessence-type fields. This 
would introduce quantitative differences, but would not change the basic qualitative 
results we are discussing. 

In the following section we describe our numerical simulations of the evolution 
of the domain wall network. We then proceed to discuss the basic features of the 
structure formation mechanism for this scenario in section Finally we present 
our results in section IV and discuss our conclusions in section |v|. 
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II. EVOLUTION OF THE DOMAIN WALLS 



We consider the evolution of a network of domain walls in a k=0 Friedmann- 
Robertson- Walker universe with line element: 

ds 2 =a 2 (r))(dr] 2 - dx 2 - dy 2 - dz 2 ), (1) 

where 0(77) is the cosmological expansion factor and 77 is the conformal time. The 
dynamics of a scalar filed <fi is determined by the Lagrangian density, 

£ = -^(^,«<r + W), (2) 

where we will take V{4>) to be a generic <fr potential with two degenerate minima 



^2 
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V(4>) = V ^0-1 +C(<f>). (3) 



where C{4>) smoothly interpolates between V_ at <fi = —(f>o and V + at 4> — +^0 and is 
otherwise such that the potential V(<f>) has two minima at ±0 O which have different 
energies. This obviously admits domain wall solutions [33| . The precise form of 
the function C(<j)) is not important as it will not affect domain walls dynamics if 
C((j}) <§C Vo for all possible values of the scalar field 0. By varying the action 

J3 



S = J dr) J d^x^/^C, (4) 
with respect to <fi we obtain the field equation of motion: 

drj 2 a drj 

with 



2 1 2 - : 3r - y2 ^ = - fl2 irr- ( 5 ) 



rf d 2 d 2 d 2 

V = 1 1 . (6) 

dx 2 dy 2 dz 2 

When making numerical simulations of the evolution of domain wall networks (or 
indeed other defects) it is also often convenient to modify the equation of motion 
for the scalar field <j) in such a way that the comoving thickness of the walls is fixed 
in comoving coordinates. This is known as the PRS algorithm |Q and it will not 
significantly affect the large-scale dynamics of domain walls. 

Hence, we will modify the evolution equation for the scalar field <j> according to 
the PRS prescription: 

ft+A-^-V^ = -^^. (7) 
or] z a or/ o<p 

where f3\ and 02 are constants. We choose 02 — in order for the walls to have con- 
stant comoving thickness and [3% = 3 by requiring that the momentum conservation 
law for how a wall slows down in an expanding universe is maintained |34[ . 

We perform two-dimensional simulations of domain wall evolution for which 
d 2 (j)/dz 2 — 0. These have the advantage of allowing a larger dynamic range and 
better resolution than tree-dimensional simulations. We assume the initial value of 
</> to be a random variable between — (f>o and cf>o and the initial value of cj> to be equal 
to zero everywhere. We normalise the numerical simulations so that cf>o = 1. We 
have checked jj5| that the initial conditions are unimportant (as expected) , because 
the domain wall network rapidly approaches a scaling solution with its statistical 
properties being independent of the initial configuration. 
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III. EVOLUTION OF COSMOLOGICAL PERTURBATIONS 



As described previously, we shall assume that the universe is made up of several 
regions (domains), with different vacuum energy densities. We will discuss the case 
where there are two such possible values, but it would be easy to generalise this to a 
distribution with a continuous range of values. Moreover, one assumes that the thin 
region separating any two of the domains considered (domain wall) is not relevant 
for structure formation. This happens if the potential of the field is small enough 
at the origin. Vacuum energy becomes dominant only for recent epochs and so we 
shall be concerned with the evolution of perturbations only in the matter-dominated 
era, neglecting the contribution of the radiation component. The average vacuum 
density is p A = (V + + V-)/2, so that 

— V + + V- 

n A = p A / Pc = — , (8) 

2,0c 

where p c is the critical density. We define A£l A in a particular domain as 

2V - W - V 

An A = A PA / Pc = -± (9) 

where V is the value of the vacuum energy density inside the domain. In our case 
Afl A can have one of two possible values Ail A = ±|V+ — V- |/2p c . 

In the synchronous gauge, the linear evolution equation for cold dark matter 
density perturbations, 5 m , in a flat universe with a non-zero cosmological constant 
can be written as 

s m + ns m - ^n 2 (n m s m - 2Ar! A ) = o, (io) 

where the evolution of the scale factor a, is governed by the Fricdmann equation 

H 2 = Hl(Tf m a- l + (l-T? m )a 2 ). (11) 

Here a dot represents a derivative with respect to conformal time, the superscript 
'0' means that the quantities are to be evaluated at the present time, and H = a/ a. 
Note that the average matter and vacuum energy densities at an arbitrary epoch 
can be written as 

- K 

O m = nn (12) 

and Q A = 1 — O m , where we have also chosen clq = 1. We start the simulation 
sufficiently early (say at a red-shift Zi = 10 3 ) so that we do not have to worry about 
the initial compensation. Because our aim is to investigate the average equation of 
state of each domain we shall assume the following initial conditions for eq. (flO|): 

& m {rii) = 0, S m {Vi) = 0. (13) 

We will parametrise the density perturbations 5 m in each domain as a fluctuation 
A£l m = pm&ml Pc = S m n m in the local value of the matter density. It is easy to 
show that the time component 

r 00 = -^H 2 (An m + A 57a) + (14) 

of the pseudo-stress-energy tensor, t^, must be compensated on super-horizon 
scales (with too oc fc 4 ) Here, the quantity —5/3TC can be interpreted as a 

fractional variation in the local expansion rate parametrised by h. We verified that: 
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Afi A + Afi ro - 2- 0, 

h 



(15) 



where Ah/h = —5/ZTL except near the boundary between different domains as 
expected from the previous discussion. Note that the factor of two in the 'Hubble' 
term of the above expression arises because the Friedmann equation relates 7i 2 to 
the average matter and vacuum energy densities. 

This can be confirmed in Fig. 1, which shows the value of the three terms above, 
as well as of its sum, for a particular simulation. 

The cancellation is almost perfect everywhere except where there are domain 
walls. A more detailed study also shows that at the domain walls the area where 
the sum of the three terms is positive is equal to that where the sum is negative, so 
that these average out to zero over the whole box. 



IV. RESULTS AND DISCUSSION 



One key assumption of our model is that the energy scale of the phase transi- 
tion which produced the domain walls is sufficiently low so that the domain walls 
have a negligible effect on structure formation. The standard bound of 1 MeV p3fl 
obviously applies here. The main consequence of this assumption is that the in- 
homogeneities are not generated by the domain walls but are due to the different 
vacuum densities in different domains. The fact that the vacuum densities only 
become important at late times explains why these inhomogeneities in the local 
values of the cosmological parameters fi A , and h are created only at recent 
times. Different values of fi A and lead to different linear growth factors from 
early times to the present. For primordial perturbations in a flat KCDM model, 
the quantity 

<?( Mm) 2 (^0+ 14Q 14Q +Vm) J W 

provides a very good fit to the suppression of growth of density perturbations rel- 
ative to that of a £l° n = 1 universe |3?],[58|]. This rescaling was also shown to be 
valid for generic topological defect models for structure formation on all scales of 
cosmological interest and for any reasonable combination of the cosmological pa- 
rameters h, and Q A |3^]. Our model will thus lead to a universe made up of 
several domains, in which the growth factor g(f^) has different values. The result 
is an inhomogeneous universe today with, for example, the abundance of clusters 
of galaxies varying from one position to another. Hence, we expect the presence 
of occasional 'great walls' separating domains with different values of the cosmo- 
logical parameters at high red-shift. However, given that in our model the domain 
walls do not generate any cosmologically relevant fluctuations these 'great walls' 
will simply provide a smooth transition between domains with different values of 
the cosmological parameters. This can clearly be seen in Figs. 1 and 2. 

Given that the length-scale corresponding to these inhomogeneities is expected to 
be close to the horizon scale (simply due to the dynamics of the domain walls), we 
may not be able to realize that we live in an inhomogeneous universe. This happens 
because as we look far away we are also looking backwards in time and the universe 
will get more and more homogeneous as the red-shift of the cosmological objects we 
are looking at increases (note that AQa = Afi^/(O m (l + z) 3 + (1 — fl m )) decreases 
very rapidly with red-shift). In fact, the higher the red-shift we are looking at 
the larger is the possibility of finding a domain with different local values of the 
cosmological parameters but the harder is to realize that. This is clearly illustrated 
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in Fig. 2, which shows the values of f2 m and Ah/h for a particular simulation. 
One should notice the different colour schemes being used for z — and 2 = 1; had 
we used the same colour scheme for the earlier redshift, no significant fluctuations 
would be detectable by visual inspection. This also explains why CMB fluctuations 
created at last scattering will be completely negligible. Significant CMB fluctuations 
can only be created at late times but even these are expected to be small if we do 
not live near the edge of our domain. However, we know that to be true (if A£!° is 
not too small) because otherwise the Universe would look very anisotropic. 

This makes the observational detection of this effect somewhat non-trivial. The 
best way of doing it should be through the determination of the number density 
of objects as a function of redshift in different directions, assuming that one has 
a reliable understanding of other possible evolutionary effects. Specific examples 
would be the counting of X-ray or Sunyaev-Zel'dovich galaxy clusters as a function 
of the red-shift z such as can be performed by XMM or Planck [ftO 41 1, large- 
scale velocity flows jl2| or gravitational lensing statistics of extragalactic surveys 
Another possibility is to look for cosmological anisotropies out to z ~ 1 with 
Supernovae la. Results from a recent analysis [Q using a combined sample of 79 
high and low red-shift supernovae are consistent with an homogeneous and isotropic 
universe but do not exclude the existence of significant anisotropies on cosmological 
scales. Unfortunately, some of their assumptions and results do not apply to our 
model. Such an investigation within the scope of our model is more complex because 
the allowed variations on the local values of the cosmological parameters are not 
independent and the spatial distribution of the patches with different values of the 
cosmological parameters is unknown. A simplified analysis of Supernova and CMB 
data constraints was recently performed in |45|,ffr| 



V. CONCLUSIONS 



In this paper we have provided a simple example of a cosmological scenario where 
the universe becomes inhomogeneous at a very recent epoch, in a way which is per- 
fectly consistent with current observations assuming that we are not very close to 
one of the boundaries. The inhomogeneity arises due to the onset of vacuum en- 
ergy domination, if the value of the 'cosmological' constant is different in different 
domains. This in turn implies that the subsequent dynamics of each patch of the 
universe will be different, leading to different values of other cosmological parame- 
ters in each patch, such as the matter density and the Hubble constant. The size of 
each patch is determined by the dynamics of the domain walls but is expected to be 
of the order of the horizon (although causality prevents it from exceeding it). This 
fact makes the detection of the cosmological signatures of this kind of model more 
difficult given that the contribution of the vacuum energy density rapidly becomes 
negligible with increasing red-shift. 
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FIG. 2. The values of Qa = + AOa (top panels) , fl m — fl m + A£l m (middle panels) 
and Ah/h (bottom panels) for a particular simulation box, for 2 = (left) and z = 1 
(right). Notice the different colour scales in the left and right panels. The horizon size at 
z — and z = 1 is respectively about 1/4 and 1/5 of the box size. 
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